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We report a search for the decays B° -> K* 1(* and B° K*°K*°. We also measure other 
charmless decay modes with K + ir~ K~ n + and K + ir~ K + ir~ final states. The results are obtained 
from a data sample containing 657 x 10 6 BB pairs collected with the Belle detector at the KEKB 
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PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



The charmless decay B° -)■ K*°K*° [l[ proceeds a K 
through electroweak and gluonic b — > d "penguin" loop 
diagrams. It provides an opportunity to probe the dy- 
namics of both weak and strong interactions, which play 
an important role in CP violation phenomena. For a B 
meson decaying to two vector particles, B — > VV, the- 
oretical models based on the frameworks of either QCD 
factorization or perturbative QCD predict the fraction 
of longitudinal polarization (/l) to be ~ 0.9 for tree- 
dominated decays and ~ 0.75 for penguin-dominated de- 
cays 0, 01]. However, the measured polarization fraction 
in the pure penguin decay B — » (j)K * has a somewhat 
lower value of /l ~ 0.5 [41. This unexpected result has 
motivated further studies [5j. 

One resolution to this puzzle is a smaller B — > K* form 
factor that could reduce /l significantly @. If this ex- 
planation is correct, the penguin-dominated decay B° — > 
K*°K*° should exhibit a similar polarization fraction. A 
time-dependent angular analysis of B° — > K*°K*° could 
distinguish between penguin annihilation and rescatter- 
ing as mechanisms for the /l observed in B — > <f>K* de- 
cays 0. The B° -)■ A'* A* mode can also be used 
to extract the branching fraction corresponding to the 
longitudinal helicity final state, determine hadronic pa- 
rameters for the b —> s decay B s — > K*°K*°, and help 
constrain the angles 4>2 (a) and ^3 (7) of the Cabibbo- 
Kobayashi-Maskawa unitarity triangle . The topologi- 
cally similar decay B° — > K*°K*° is strongly suppressed 
in the Standard Model (SM) ; its observation would indi- 
cate new physics. 

Theoretical calculations predict the branching frac- 
tions for B° -> K*°K*° and B° K*°K*° to be (0.17- 
0.92) x 10~ 6 @ and (2.9 ±0.2) x 10" 15 0, respectively. 
The BaBar collaboration has reported a branching 
fraction B{B a K*°K*°) = (1.28+°;^ ± 0.11) x 10~ 6 
and has set an upper limit B(B° -» K*°K*°) < 0.41 x 
10~ 6 at the 90% confidence level (C.L.). In this pa- 
per, we report a search for the decays B° — > K*°K*°, 
B° — > K*°K*°, and other charmless decay modes with B c 



or K+tt-K^ 



final state using a data 



sample 1.7 times larger than that of BaBar. The data 
sample used in the analysis contains 657 x_ 10 6 BB 
pairs collected with the Belle detector 



12, hj at the 



KEKB asymmetric-energy e + e (3.5 GeV and 8 GeV) 



collider [14( operating at the T(AS) resonance. The Belle 
detector includes a silicon vertex detector, a 50-layer cen- 
tral drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), and a barrel-like arrange- 
ment of time-of- flight scintillation counters (TOF). Sig- 
nal Monte Carlo (MC) events are generated with EVT- 
GEN [l5[ , and final-state radiation is taken into account 
with the PHOTOS package Generated events arc 

processed through a full detector simulation program 
based on GEANT3 " 
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We reconstruct signal decays from neutral combina- 
tions of four charged tracks fitted to a common vertex. 
Neutral K* mesons are reconstructed via K*° — »• K + ir~ 
and A"* — > K~tt + . Charged track candidates are re- 
quired to have a distance of closest approach to the in- 
teraction point of less than 2.0 cm in the direction along 
the positron beam (z axis), and less than 0.1 cm in the 
transverse plane. Tracks are also required to have a lab- 
oratory momentum in the range [0.5, 4.0] GeV/c, a polar 
angle in the range [32.2, 127.2]°, and a transverse momen- 
tum pt > 0.1 GeV/c. Charged pions are identified using 
information from the CDC (dE/dx), ACC, and TOF de- 
tectors 



18|. 



We distinguish charged kaons from pions 
using a likelihood ratio TZk = £k/(£-k + where 
£-n{C-K) is the likelihood value for the pion (kaon) hy- 
pothesis. We require TZk < 0.4 (TZk > 0.6) for the two 
charged pions (kaons). The kaon (pion) identification ef- 
ficiency is 83% (90%), and 6% (12%) of pions (kaons) 
are misidentified as kaons (pions). We also use the lep- 
ton identification likelihood TZ X (x denotes cither fi or 
e) described in Ref. jl^]: charged particles identified as 
electrons (TZ e > 0.9) or muons (TZ^ > 0.9) are removed. 
We veto B -> D^X, B -» D±X, B D°X, and 



4>X decays that result in the A + 7r K 



final 
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state, and we veto B -> D^±X and B -> £>°A de- 
cays that result in the if + 7r~i ; r + 7r~ final state. For the 
B -> D^^X and B -> vetos, we remove can- 

didates that satisfy either \M{K ± K^tt t ) — m flT | < 

13 McV/c 2 , |M(K ± 7rT 7 rT) - m DT | < 13 McV/c 2 , 

|M(lf ± ^7r : ' : )-m D = F | < 13MeV/c 2 , or \M{K^h^)- 
m D ^ \ < 13 MeV/c 2 , where m fl T are the masses of the 

mesons, and h J- (/ij) is the kaon (pion) mass as- 
signed to a pion (kaon) candidate track [i.e., a kaon 
(pion) was misidentified as a pion (kaon)]. For the 
B —> D a X veto, we remove candidates satisfying ei- 
ther \M(K ± K T ) - m D o \ < 13 McV/c 2 or \M{K ± h^) - 
m D o\ < 13 MeV/c 2 , where m D o is the mass of the D° 
meson. For the B — > <pX veto, we remove candidates 
satisfying \M(K ± h] { ) - m^,| < 20 MeV/c 2 , where is 
the mass of the <f> meson. These vetos together remove 
9.7% (3.6%) of longitudinally polarized B° -> K*°K*° 
(B° — > K*°K*°) signal, according to MC simulation. 

Signal event candidates are characterized by two 
kinematic variables: the beam-energy-constraincd mass, 
M hc = v^bcam ~~ P B' and tlie energy difference, AE = 
E* B — i^bcam, where E'beam is the run-dependent beam en- 
ergy, and Pg and Eg are the momentum and energy of 
the B candidate in the T(45) center-of-mass (CM) frame. 
We distinguish nonresonant B° —> KKirir decays from 
our signal modes by fitting the two-dimensional mass 
distributions M(X + 7r _ ) vs. M(K~n+) or M(K+tt~) 
vs. M(K + n~). There are two possible combina- 
tions in B° -4 K*°K*° reconstruction for M(K + ir~) 
vs. M(K+n-): (K+^)(K+^) and {K+^)(K+*p, 
where the subscripts label the momentum ordering, i.e., 
Ky (tTi) has higher momentum than K^n^)- We con- 
sider both (K^tt^)(K2^2) an d (^i" 7r 2~)(^2 h7r r) combi- 
nations and select candidate events if either one of the 
combined masses lies within the signal window of [0.7, 
1.7] GeV/c 2 . If both combinations fall within the signal 
window, we select the (K^tt^)^^^) combination. Ac- 
cording to MC simulation, this choice selects the correct 
combination for signal decays 99% of the time. For fit- 
ting, we symmetrize the M 2 (K + tt~) vs. M 2 (K + tt~) plot 
by plotting M 2 (K^tt^) [M 2 (K^tt~)} on the horizontal 
axis for events with an even [odd] event number. This 
number denotes the location of the event in the data set 

(i.e., n eve nt = l,2,3..JV t otal)- 

The dominant source of background is continuum 
e + e _ — > qq (q = u,d,s and c) events. To distinguish 
signal from the jet-like continuum background, we use 
modified Fox- Wolfram moments that are combined 
into a Fisher discriminant. This discriminant is subse- 
quently combined with the probabilities for the cosine of 
the B flight direction in the CM frame and the distance 
along the z axis between the two B meson decay vertices 
to form a likelihood ratio 1Z = £ s /(£ s + C q q). Here, 
C s {Cqq) is a likelihood function for signal (continuum) 



events that is obtained from the signal MC simulation 
(events in the sideband region M bc < 5.26 GeV/c 2 ). For 
additional suppression, we also use a flavor tagging qual- 
ity variable r provided by the Belle tagging algorithm [2(J 
that identifies the flavor of the accompanying B° meson 
in the T(45) —5- B°B° decay. The variable r ranges from 
r = for no flavor discrimination to r = 1 for unambigu- 
ous flavor assignment, and it is used to divide the data 
sample into six r bins. As the discrimination between 
signal and continuum events depends on the r-bin, we 
impose different requirements on 1Z for each bin. The 
requirements are determined by maximizing a figure-of- 
merit N s / y/W a + N qq , where N s (N qq ) is the expected 
number of signal (continuum) events in the signal region 
AE g [-0.045, 0.045] GeV, M bc g [5.27, 5.29] GeV/c 2 , 
and Mi, 2 (l£V) g [0.826, 0.966] GeV/c 2 . 

In about 17% of events there are multiple B° — > 
K*°K*° or B° -> K*°K*° candidates. For these events 
we select the candidate with the smallest \ 2 value for 
the B° decay vertex reconstruction. This selects the 
correct combination 87% (97%) of the time for longi- 
tudinally (transversely) polarized B — > K*°K and 
B° — > K*°K *° decays. The overall reconstruction ef- 
ficiency for B° — > K*°K*° as obtained from MC sim- 
ulation is 4.43% (5.23%) for longitudinal (transverse) 
polarization. The overall reconstruction efficiency for 
B° -> K*°K*° is 5.74% (5.92%) for longitudinal (trans- 
verse) polarization. The efficiency for longitudinal polar- 
ization is lower, as in this case the K*° daughters produce 
lower momentum kaons and pions. 

The signal yields for B° -> K*°~K*° and other B° -)• 
K + tt~ K~tt + decays are extracted by performing an ex- 
tended unbinned maximum likelihood (ML) fit to the 
variables M bc! AE, M x , and M 2 , where Mi = M(K + n~) 
and M 2 = M(K~tt + ). This four-dimensional fit dis- 
criminatesamong K*°K*°, K*°Kir, i^(1430)iTg(1430), 
K${U30)K*°, K£(U30)Ktt, and nonresonant KKttw fi- 
nal states. Since there are large overlaps between these 
states, we distinguish them by fitting a large (Mi,M 2 ) 
region: M ij2 G [0.7, 1.7] GeV/c 2 . We use a likelihood 
function 

£ = exp(-2>) ff(5>jA (1) 

3 »=1 V 3 J 

where i is the event identifier, j indicates one of the 
event type categories for signals and backgrounds, rij 
denotes the yield of category j, and Pj is the prob- 
ability density function (PDF) of event i for cate- 
gory j. The PDF is a product of two smoothed 
two-dimensional functions: Pj = Pj(Ml c ,AE l ) x 
r,i .\l! ..\iy : ee PjiMl^AE^MiMZ). The signal yields 
for B° -> K*°K*° and other B° -> K+tt-R+tt- decays 
are extracted by another four-dimensional fit in the same 
way, except that for this fit M 2 ee M(K + ir~). 

For the B signal components, the smoothed functions 
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P(M bc , AE) and P(M 1 ,M 2 ) are obtained from MC sim- 
ulation. For the Mb c and AE PDFs, possible differences 
between data and the MC modeling are calibrated using 
a large control sample of B° — > D~ (K + tt~tt~)tt + de- 
cays. The signal mode PDF is divided into two parts: 
one is correctly reconstructed events (CR), and the other 
is "self-cross-feed" events (SCF) in which at least one 
track from the signal decay is replaced by one from the 
accompanying B decay. We use different PDFs for CR 
and SCF events and fix the SCF fraction (/scf) to that 
obtained from MC simulation, i.e., 

Signal = (1 - /SCF) X PcrMc, MlM) 

+ ZscF x P S cv(ML AE\M{, Ml) . (2) 

For the continuum and b — > c decay backgrounds, we 
use the product of a linear function for AE, an ARGUS 
function 21 1 for Mb c , and a two-dimensional smoothed 
function for M 1 -M 2 - The shape parameters of the linear 
and ARGUS functions for the continuum (b — > c) events 
are floated (fixed) in the fit; the shape of the Mi~M 2 func- 
tions for the continuum and b — > c events are obtained 
from MC simulation and fixed in the fit. The yields of 
the continuum and b — > c decay backgrounds are floated 
in the fit. For the charmless B decay backgrounds, we 
use separate PDFs for B° — >• K + K~ K° , nonresonant 
B° — > Ktttttt, and other charmless B modes; all the PDFs 
are obtained from MC simulation. Note that the nonres- 
onant B a — > Ktttttt decay will enter the sample if one 
of the pions is misidentified as a kaon; in this case the 
mean of the AE distribution shifts by about +75 McV, 
since assigning a kaon mass instead of a pion mass in- 
creases the B candidate energy. In the fit, we fix the 
yield of B° — > K + K~K° to 32 events, corresponding to 
a branching fraction of 24.7 x 10 -6 [22j, and the yield of 
other known charmless B decays to that expected based 
on world average branching fractions j23|. We set the 
branching fraction for B° — > K 2 (1430)X to zero and only 
consider it in the systcmatics, as this mode has a large 
correlation with B° — > Kq( 1430) X. The yield of nonres- 
onant B° — > Ktttttt is floated. For the fully nonresonant 
modes, we assume the final-state particles arc distributed 
uniformly in three- and four-body phase space. 

The fit results are listed in Table HI and projec- 
tions of the fit superimposed to the data are shown in 
Figs. Q] and [2j The statistical significance is calculated 
as yj — 2 ln(£o/ £max ) , where Co and £ max are the val- 
ues of the likelihood function when the signal yield is 
fixed to zero and when it is allowed to vary, respectively. 
We do not find significant signals for B° -t K*°K*°, 
B° -> 

K+TT- 



K* U K* U , and other charmless decay modes with 
K^tt final states, and determine 90% C.L. upper 
limits for the yields (N). These limits are calculated via 



Jo_ C(x)dx 
Jo°° C(x)dx 



where x corresponds to the number of signal events. We 
include the systematic uncertainty in the upper limit 
(UL) by smearing the statistical likelihood function by 
a bifurcated Gaussian whose width is equal to the to- 
tal systematic error. We also smear C when calculating 
the signal significance, except that only the additive sys- 
tematic errors related to signal yield are included in the 
convolved Gaussian width. Our upper limits correspond 
to a longitudinal polarization fraction /l = 1; as the ef- 
ficiency for /l < 1 is higher than that for /l = 1, our 
limits are conservative. 

To check our reconstruction efficiencies, we mea- 
sure the yields of control samples B° — > D~tt + — > 



(K+K~ 



r+ and B° -> D°K* 



(K + TT-)(K+TT-) 



These modes have a similar topology to the signal modes 
and are selected using the same selection criteria ex- 
cept that, instead of D vetos, we require \M(KKtt) — 
m D ± | < 13 MeV/c 2 , \M (Ktt) - m D o | < 13 MeV/c 2 , and 
826 MeV/c 2 < M(Ktt) < 966 MeV/c 2 . The efficiencies 



are 19% for B c 



D-tt+ and 11% for B c 



'0 if*0 



D°K 



0.90. 



(3) 



The difference in signal yields between the measured and 
expected values are (5.8 ±5.8)% and (5.6 ± 27.8)% for 
B" -> D-tt + and B° -> D°K*°, respectively. These 
differences are consistent with zero. 

The systematic errors (in units of events) are sum- 
marized in Tables [TT] and Mil For systematic uncer- 
tainties due to fixed yields, e.g., that of charmless B 
background, we vary the yields by their uncertainties 
(±lcr). For the systematic uncertainties due to B° —> 
isT 2 *(1430)X decays, including B° — » ^1430)^3(1430), 
B° -> if*(1430)^(1430), B° -> ^(1430)if*°, and 
B° -> K%(1430)Ktt, we float their yields in the four- 
dimensional ML fit; the differences between these results 
and the nominal fit values are taken as systematic er- 
rors. Systematic uncertainties for the AE-M^ C PDFs 
arc estimated by varying the signal peak positions and 
resolutions by ±lcr and repeating the fits. Systematic 
uncertainties for the M\-M 2 PDFs are estimated in a 
similar way; we vary the mean and width of K*° and 
fTo(1430) mass shapes according to the uncertainties in 
the world average values 22| . A systematic error for the 
longitudinal polarization fraction is obtained by chang- 
ing the fraction from the nominal value /l = 1 to the 
lowest possible value /l = when evaluating the recon- 
struction efficiency. According to MC simulation, the 
signal SCF fractions are 13.4% for (longitudinally polar- 
ized) B° -> K*°K*°, 7.9% for B° -> K*°Ktt, 6.7%_for 
B° -)• ^(1430)7^(1430), 6.7% for B° -> Jf *( 1430)7?*°, 
7.6% for B° ->■ KS(U30)Ktt, and 9.2% for nonresonant 
B° — > KKtttt. We estimate a systematic uncertainty due 
to these fractions by varying them by ±50%. 

A high-statistics MC study indicates that there are 
small fit biases; these are listed in Table |U We find that 
fit biases occur due to the correlations between the two 
sets of variables (AE, iWb c ) and (Mi, M 2 ), which are 
not taken into account in our fit. We correct the fitted 
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TABLE I: Fit results for decay modes with final states K + tt~ K~iv + and K + ir~ K + ir~ . The fit bias (in units of events) is 
obtained from MC simulation; the yield includes the bias correction; the efficiency e includes the PID efficiency correction and 
branching fractions for A* — > K + tv~ and Aq(1430) — ► K + ir~ (66.5% and 66.7%, respectively); and the significance 5 is in 
units of a. The first (second) error listed is statistical (systematic). 
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< 0.8 


B° -> K*°K-n+ 




-5.4 ±2.9 
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FIG. 1: Projections of the four-dimensional fit onto (a) AE, (b) Mb c , (c) M(K + n ), and (d) M(K n + ) for candidates 
satisfying (except for the variable plotted) AA G [-0.045, 0.045] GeV, M bc G [5.27, 5.29] GeV/c 2 , and AA,2(Att)_G 
[0.826, 0.966] GeV/c 2 . The thick solid curve shows the overall fit result; the solid shaded region represents the B° — > K*°K*° 
signal component; and the dotted, dot-dashed and dashed curves represent continuum background, b — > c background, and 
charmless B decay background, respectively. 




AE (GeV) M bc (GeV/c 2 ) 



FIG. 2: Same as for Fig. □ but for the B° -> K*°K*° (A 
M 2 (A+7r-). 




°J.7 0.9 1.1 1.3 1.5 1.7 Vl.7 0.9 1.1 1.3 1.5 1.7 

M,(KV) (GeV) M 2 (KV) (GeV) 



7r")(A+7r") study: (a) AA, (b) M bc , (c) Mi(K + n~), and (d) 
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yields for these biases. To take into account possible dif- 
ferences between MC simulation and data, we take both 
the magnitude of the bias corrections and the uncertainty 
in the corrections as systematic errors. The systematic 
errors for the efficiency arise from the tracking efficiency, 
PID, and the 1Z requirement. The systematic error on 
the track- finding efficiency is estimated to be 1.2% per 
track using partially reconstructed D* events. The sys- 
tematic error due to the PID is 1.0% per track as esti- 
mated using an inclusive D* control sample. The sys- 
tematic error for the 1Z requirement is determined from 
the efficiency difference between data and MC samples of 
B° ->■ D~(K + TT~n~)n + decays. 

In summary, we have used a data sample correspond- 
ing to 657 x 10 6 BB pairs to search for B° K*°K*°, 
B° —> K*°K*°, and other charmless decay modes with 
a K + tt~ A T 7r ± final state. We do not find significant 
signals for any of these modes. Our measured branching 
fraction for B° -> K*°K*° is (0.26^;^^;^) x 10" 6 , 
which is lower than that obtained by BaBar [ll| by 
2.2cr. Our 90% C.L. upper limits arc 0.8 x 10" 6 for 
B(B° -> K*°K* a ) and 0.2 x 10~ 6 for B{B° -> K*°K*°); 
those for other decay modes are listed in Table HI 
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TABLE II: Summary of systematic errors (in units of events) for decay modes with a final state K + ir~ K~ir + . The parameters 
^s°^if*(i430)x an d Nb- iUit i,s (the other known charmless B decays) correspond to branching fraction uncertainties for these 
charmless B decays. Values for /l and /scf are the uncertainties for longitudinal polarization and self-cross-feed, respectively. 



Source 


A* U A*° 


A* u A",t+ 


A * (1430) Kl (1430) A * (1430) A* u 


A *(1430)A'- 


tt + Nonresonant K + tt K n + 


Fitting PDF 


±1.8 


±40.3 


±55.6 


±37.7 


±158.0 


±102.4 


-^B0->X|(1430)X 


+1.1 


±10.2 


-7.1 


-20.4 


-52.8 


-10.5 




±0.0 


±0.1 


±0.2 


±0.1 


±0.7 


±1.0 


h 


-0.1 












/scf 


±0.7 


±1.4 


±5.6 


±2.2 


±17.2 


±14.6 


Fit Bias 


+ 1.5 
-0.7 


+2.9 
-5.4 


±5.1 


+ 13.3 
-2.3 


+8.9 
-14.7 


+7.3 
-10.8 


Tracking 


±0.4 


±0.8 


±3.5 


±0.9 


±9.6 


±6.8 


PID 


±0.4 


±0.7 


±2.9 


±0.7 


±8.2 


±6.0 


1Z requirement 


±0.2 


±0.4 


±1.6 


±0.4 


±4.5 


±3.2 




±0.1 


±0.3 


±1.1 


±0.3 


±3.1 


±2.2 


Sum 


+2.8 
-2.1 


+41.7 

-40.7 


+56.3 

-56.8 


+40.0 

-43.0 


+ 159.7 

-168.6 


+ 104.1 

-104.9 



TABLE III: Same as for Table ITT1 but for decay modes with a final state A + tt A" + tt 



Source 


A*°A* U 


A*°A+7r- 


A *(1430)Ao* 


(1430) A *(1430)A* U 


Nonresonant A + 7r A + 7r 


Fitting PDF 


±2.4 


±40.0 


±20.7 


±22.5 


±30.5 


Nb° ->K| (1430)X,K^ (1430)Xit 


-0.3 


±16.8 


±85.1 


-20.2 


-96.8 




±0.0 


±0.4 


±0.1 


±0.3 


±0.3 


/l 


-0.7 










/scf 


±0.1 


±0.2 


±0.8 


±0.5 


±0.2 


Fit Bias 


+0.5 
-1.0 


±2.7 


+ 1.3 
-3.4 


+8.2 
-1.6 


+ 7.7 
-2.2 


Tracking 


±0.2 


±0.0 


±1.2 


±0.4 


±0.5 


PID 


±0.2 


±0.0 


±1.0 


±0.3 


±0.4 


TZ requirement 


±0.1 


±0.0 


±0.6 


±0.2 


±0.2 


N B b 


±0.1 


±0.0 


±0.4 


±0.1 


±0.2 


Sum 


+2.5 
-2.7 


+43.5 

-40.1 


+87.7 
-21.1 


+ 23.9 

-30.3 


+31.4 

-101.5 



